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1.  Brief  Progress  Summary  and  Discussions  of  Published  Manuscripts. 

We  have  investigated  properties  of  III-V  compound  lattice  matched  heterolayers  to  explore 
their  ultimate  limitations  and  advantages  for  high  speed  semiconductor  devices.  During  the  past 
three  years  of  this  grant,  twelve  manuscripts  (abstracts  attached)  were  published  or  accepted  for 
publication  and  eight  invited  presentations  have  been  given.  The  three  most  significant  achieve¬ 
ments  have  been  the  development  of : 

(i)  The  first  Monte  Carlo  simulation  of  the  high  electron  mobility  transistors  (HEMT) 

(ii)  The  first  Monte  Carlo  model  of  planar  doped  barrier  PDB  transport  including  plasmon 
scattering  and  the  Pauli  principle,  and 

(iii)  A  new  ultrafast  switching  device  based  on  a  combination  of  tunneling  and  hot  electron 
thermionic  emission. 

Manuscript  [l]  summarizes  high  speed  transport  in  submicron  high  electron  mobility  transis¬ 
tors  (HEMT's).  The  manuscript  deals  mainly  with  the  physical  effects  which  are  important  for  the 
HEMT  operation  and  to  a  lesser  extent  with  actual  device  performance. 

The  simulations  show  that  significant  velocity  over-shoot  exists  in  these  structures  and  precise 
values  of  the  over-shoot  velocity  are  given.  The  device  speed  is  substantially  enhanced  by  this 
effect  for  devices  of  1/2  micrometer  feature  sizes  and  below. 

We  have  completed  the  simulation  code  for  PDB's  including  standard  transport  effects  as  well 
as  scattering  by  coupled  phonon-plasmons.  electron-electron  pair  interactions,  ionized  impurity 
scattering  and  also  the  Pauli  principle. 

The  results  of  our  calculations  have  been  compared  with  experiments  of  Hayes  and  coworkers 
at  AT&T  Bell  Laboratories  and  of  Heiblum  and  coworkers  at  IBM. 

We  found  that  there  is  an  overall  agreement  between  theory  and  experiment  (manuscript  II) 
but  the  experimentally  found  energetic  width  (broadening)  of  the  ballistic  peak  is  much  greater 
than  predicted  by  the  theory  in  the  Hayes  experiment  and  slightly  greater  in  the  experiment  by 
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Heiblum  et  al.  We  attribute  the  broadening  of  the  Hayes  experiment  to  the  fluctuating  impurity 
potentials  at  the  barriers  (there  are  only  a  few  impurities  across  the  width  of  thin  barriers). 
GaAs  —Alx  Ga  i_y  As  heterolayers  and  therefore  their  barriers  do  not  suffer  from  these  fluctuations. 
The  broadening  in  the  experiments  of  Heiblum  et  al.  can  be  attributed  to  electron-electron  scatter¬ 
ing.  A  report  about  the  experimental  work  and  our  theory  has  appeared  in  Science  (Vol.  231. 
p.22).  Details  are  given  in  the  attached  manuscript. 

Manuscript  (III)  demonstrates  a  new  semiconductor  switching  phenomenon  and  infers  from  a 
phenomenological  theory,  switching  times  of  ^  200  femtoseconds.  We  are  especially  excited  about 
this  project  which  could  lead  to  an  enormously  fast  semiconductor  switching  circuit  element. 

Perhaps  the  most  complete  study  of  interface  electronic  transport  in  (5)  subbands  at  the 
AlGaAs-GaAs  interface  has  been  started  by  Yokoyama  (manuscript  IV.  V)  and  is  currently  being 
completed  by  Artaki  who  also  includes  transport  at  low  temperatures  (4.2K). 

Smaller  projects  have  included  time  dependent  ensemble  Monte  Carlo  simulations  for  planar 
doped  GaAs  structures  (manuscript  VI)  transport  in  superlattices  including  collision  broadening 
(manuscript  VII)  and  effects  of  size  quantization  on  electrons  and  hole  transport  (Manuscript  VIII 
and  IX).  While  the  latter  papers  are  mainly  devoted  to  the  inclusion  of  quantum  effects  into  other¬ 
wise  semiclassical  Monte  Carlo  simulations,  we  have  also  taken  the  first  steps  to  a  full  quantum 
transport  theory  involving  Feynman  path  integrals.  The  work  which  has  been  accomplished 
includes:  Solutions  of  the  linearly  coupled  electron-acoustic  phonon  system  and  comparison  of 
results  with  analytic  work  and  calculation  of  the  density  matrix  and  expectation  values  of  elec¬ 
tronic  properties  as  a  function  of  time. 

Our  experience  in  the  projects  outlined  above  has  been  summarized  in  two  book  chapters 
(manuscript  X  and  XI)  as  well  as  in  eight  invited  presentation  (manuscript  XII  and  list  below). 

2.  Invited  Presentations: 


Principles  of  Hot  Electron  Thermionic  Emission  (Real  Space  Transfer)  in  Semiconductor 


Heterolayers  and  Device  Applications".  1st  Intern.  Workshop  on  Future  Electron  Devices. 
Tokyo.  February  1984. 

"Transient  Electron  Transport  in  Semiconductor  Heterolayers".  Gordon  Conference  on  MIS 
Systems.  Tilton.  NH.  July  1984. 

"Principles  of  Hot  Electron  Thermionic  Emission  (Real  Space  Transfer)  in  Semiconductor 
Heterolayers  and  Device  Applications."  European  Physical  Society  5th  General  Conference 
of  the  Condensed  Matter  Division  of  the  EPS.  Berlin  (West).  March  18-22.  1985. 

"Electronic  Transport  in  Semiconductor  Heterolayers  and  Promising  Device  Concepts".  Pro¬ 
gram  of  the  34th  WEH-Seminar.  Bad  Honnef.  May  13-15.  1985. 

"Properties  and  Applications  of  Multiquantum  Well  and  Superlattice  Structures."  Ameri¬ 
can  Vacuum  Society  32nd  National  Symposium,  Houston.  TX,  November  19-21.  1985. 

"Analysis  and  Simulation  of  Semiconductor  Devices."  Scientific  Applications  and  Algorithm 
Design  for  High  Speed  Computing  Workshop,  Urbana.  IL.  April  7-10.  1986. 

"Computations  of  Electronic  Transport  in  Semiconductor  Heterolayer  Structures."  The 
First  US-Japan  Joint  Workshop  on  Electronic/  Optoelectronics  with  Focus  on  "Optical  Bis¬ 
tability.  Optical  Computation  and  Quantum  Well  Devices",  Tokyo.  Japan.  May  28-30. 
1986. 

"Monte  Carlo  Simulations  of  Electronic  Transport  in  Semiconductor  Heterolayers."  The 
Conference  on  the  Forefronts  of  Large-Scale  Computational  Problems."  Troy.  NY  16-18. 
1986 

Scientific  Personnel: 

k  Hess  Professor 

J.  P.  I.eburton  Assistant  Professor 
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B.  Mason  Post  Doctoral  Research  Associate 

M.  Artaki  (Ph.D.  student:  Ph.D.  final  Oct.  1986) 

T.  Wang  (Ph.D.  student;  Ph.D.  completed  1985) 

K.  Brennan  (Ph.D.  student:  Ph.D.  completed  1984) 

4.  Interactions. 

Much  of  the  work  has  been  performed  in  close  cooperation  with  Dr.  G.  J.  lafrate.  Ft.  Mon¬ 
mouth. 
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Calculation  of  the  electron  velocity  distribution  in  high  electron  mobility 
transistors  using  an  ensemble  Monte  Carlo  method 

T.  Wang  and  K.  Hess 

Department  of  Electrical  Engineering,  and  the  Coordinated  Science  Laboratory ,  University  of  Illinois, 

Urbana,  Illinois  61801 

(Received  23  January  1983;  accepted  for  publication  6  February  1983) 

The  electron  velocity  distribution  is  calculated  for  an  idealized  model  of  the  high  electron  mobility 
transistor  using  a  many-particle  Monte  Carlo  model  and  a  self-consistent  two-dimensional 
Poisson  solver.  Hot  electron  effects,  nonstationary  effects,  and  real  space  transfer  are  analyzed. 

The  results  show  that  significant  velocity  overshoot,  2.8x  107  cm/s  at  300  K.  and  3.7  x  107  cm/s  at 
77  K  exists  under  the  gate  and  that  the  velocity  overshoot  is  limited  by  both  fc-space  transfer  and 
real-space  transfer.  The  values  of  the  overshoot  velocities  are  much  smaller  than  those  obtained 
from  the  more  conventional  drift-diffusion  model. 


I.  INTRODUCTION 


II.  METHOD  OF  ANALYSIS 


In  recent  years,  the  high  electron  mobility  transistor 
(HEMT)  has  attracted  considerable  interest  because  of  its 
potential  for  ultrahigh-speed,  low-power  logic  applica¬ 
tions.  1-4  The  key  feature  of  the  HEMT  leading  to  its  superior 
properties  and  distinguishing  it  from  other  field-effect  tran¬ 
sistors  is  the  greatly  enhanced  electron  velocity  in  the  unin¬ 
tentionally  doped  GaAs  channel.  Currently,  considerable 
theoretical  and  experimental  effort  is  underway  to  deter¬ 
mine  the  precise  reasons  for  the  speed  advantage  under  high- 
field  and  transient  conditions.5'®  It  has  been  demonstrated 
that  the  high  switching  speed  of  the  HEMT  is  due  to  the 
following  factors:  (i)  the  high  mobility  in  the  source  access 
region,  (ii)  overall  quantum  effects  in  the  conducting  chan¬ 
nel,  and  (iii)  significant  nonstationary  effects  in  the  high-field 
range  (velocity  overshoot  and  ballistic  transport).’1"111 
Among  these  factors,  the  quantum  effects  have  been  recently 
investigated  by  Widiger  etal.*  9  and  found  to  be  of  little  con¬ 
sequence  for  electrons  that  are  strongly  heated  by  the  drain 
field  and  hence  cannot  easily  be  confined  to  the  interface 
potential  well.  The  drift-diffusion  model  of  Widiger  et  al.*  is 
advantageous  for  describing  the  overall  features  of  the 
HEMT.  However,  a  more  sophisticated  approach  based  on 
the  Boltzmann  transport  equation  is  needed  to  characterize 
the  nonstationary  phenomena  which  arise  in  the  high-field 
region  close  to  the  drain  contact.  It  is  the  purpose  of  this 
paper  to  present  a  many-particle  Monte  Carlo  simulation 
including  a  complete  band  structure  for  both  the  GaAs  and 
the  AlGaAs  layers  as  well  as  all  relevant  scattering  mecha¬ 
nisms  and  a  self-consistent  two-dimensional  solution  of 
Poisson’s  equation.  Our  code  appears  to  be  well  suited  to 
explore  the  fundamental  physics  of  hot  electron  transport  in 
HEMTs,  especially  the  effects  of  k  space  and  real-space 
transfer"  and  velocity  overshoot  We  also  study  the  effect  of 
interekctronic  collisions  on  the  device  performance. 

Furthermore,  a  comparison  between  the  drift-diffusion 
model  and  the  Monte  Carlo  simulation  is  presented  and  the 
validity  of  the  drift-diffusion  equation  in  the  high  field  region 
is  examined  by  directly  comparing  overshoot  velocities  ob¬ 
tained  from  these  two  models. 


A.  Device  geometry  and  boundary  conditions 

The  device  structure  illustrated  in  Fig.  1  is  used  in  all 
the  simulations.  The  approximate  vertical  source  and  drain 
contacts  used  in  the  simulated  structure  generate  virtually 
no  difference  compared  to  the  planar  ones  in  real  devices 
since  the  potential  drop  in  the  regions  directly  under  the 
source  and  drain  is  insignificant.12  We  have  chosen  a  rela¬ 
tively  low  impurity  concentration  of  2x  1017  cm" J  in  the 
AlGaAs  layer  and  in  the  n  *  contact  regions  to  avoid  compli¬ 
cations  arising  from  the  Pauli  exclusion  principle.  Although 
such  a  low  doping  level  is  unrealistic  in  real  devices,  it  is  of 
little  consequence  to  the  main  subject  of  our  study,  the  trans¬ 
port  behavior  of  electrons  in  the  high-field  channel,  because 
the  carrier  concentration  in  the  high  field  region  is  always 
low. 

Charge  neutrality  is  assumed  in  the  regions  adjacent  to 
the  source  and  drain  contacts  throughout  the  entire  simula¬ 
tion  and  serves  as  a  boundary  condition  for  the  electron  con¬ 
centration.  The  boundary  conditions  for  the  electric  poten¬ 
tial  which  are  used  in  the  model  are  the  Dirichlet  condition 
(constant  potential)  at  all  metal-semiconductor  surfaces  and 
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FIG.  1.  Schematic  representation  of  the  computer  simulation  geometry  of 
the  AlGaAa-GaAi  HEMT.  The  impurity  concentration  is  2  X  10”  cm  ~  ‘  in 
the  n-AlOaAs  layer  and  in  n  'GaAs  contact  regions.  The  GaAs  channel  is 
unintentionally  doped.  An  effective  gate  bias  (including  the  Schottky  bar¬ 
rier  height)  of  —  0.35  V  is  applied  and  the  drain  bias  is  0.3  V. 
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Monte  Carlo  simulations  of  hot-electron  spectroscopy  in  planar-doped 
barrier  transistors 

T.  Wang  and  K.  Hess 

Department  of  Electrical  Engineering  and  the  Coordinated  Science  Laboratory.  University  of  Illinois,  Urbana, 

Illinois  61801 

G.J.  lafrate 

U.  S  Army  Electronics  Technology  and  Devices  Laboratory,  Fort  Monmouth.  New  Jersey  07703 

(Received  4  November  1985;  accepted  for  publication  27  November  1985) 

We  have  developed  a  microscopic  model  to  study  the  hot-electron-spectroscopy  method  which  * 
uses  GaAs  planar-doped  barrier  transistors.  Our  simulation  is  based  on  the  Monte  Carlo  method 
and  includes  the  effects  of  ionized  impurity  scattering,  pair  electron-electron  scattering,  long- 
range  plasmon  scattering,  and  coupled  plasmon/phonon  scattering.  The  nonparabolicity  of  the 
band  structure  and  the  Pauli  exclusion  principle  are  also  taken  into  account  in  the  highly  doped 
base  region.  The  numerical  results  show  that  the  experimental  method  of  Hayes  will  indeed  reflect 
the  overall  momentum  distribution  of  injected  hot  electrons  if  the  planar-doped  barriers  are 
“ideal.”  Ideal  means  that  the  self-consistent  potential  (due  to  conduction  electrons  and  ionized 
impurities)  is  well  described  by  the  continuum  approximation.  We  demonstrate  that  potential 
fluctuations  arising  from  the  discrete  natue  of  the  charges  and  reflection  of  electrons  at  the  base- 
collector  junction  make  it  impossible  to  obtain  the  precise  distribution  function  from  the 
experiments. 


In  recent  years,  extensive  theoretical  analyses  have  been 
performed1-4  to  assess  the  magnitude  of  velocity-overshoot 
effects  and  the  possibility  of  the  concept  of  ballistic  trans¬ 
port.3  Experimental  verifications  of  the  numerous  theories 
have  been  rare  and  very  indirect.  The  results  of  the  Monte 
Carlo  simulations  of  planar-doped  barrier  ( PDB )  devices  by 
Hesto  et  al.b  and  Wang  et  al }  suggested  the  usefulness  of 
PDBs  for  the  investigation  of  ballistic  transport  phenom¬ 
ena.7  In  a  series  of  papers,  Hayes  et  al.*~'°  presented  experi¬ 
mental  results  for  the  current  voltage  characteristics  of 
planar-doped  barrier  transistors.  They  inferred  that  the  elec¬ 
tron-momentum  distribution  can  be  obtained  by  differenti¬ 
ating  their  measured  current  with  respect  to  the  collector 
voltage.  Using  this  relation,  they  have  been  able  to  deduce 
important  results  for  the  electron  mean-free  time  of  flight, 
the  velocity  overshoot,  and  the  fraction  of  ballistic  ( unscat¬ 
tered)  electrons  which  travel  through  the  transistor  base. 

The  purpose  of  this  paper  is  twofold.  First,  we  develop  a 
microscopic  theory  of  transport  in  planar-doped  barrier  de¬ 
vices  by  using  an  ensemble  Monte  Carlo  model  which  in¬ 
cludes  all  known  scattering  mechanisms  ( including  the  elec¬ 
tron  interaction  with  coupled  plasmon/phonon  modes)  as 
well  as  important  quantum  principles  such  as  the  Pauli  prin¬ 
ciple,  which  is  not  usually  included  in  these  simulations.  Sec¬ 
ond,  we  apply  this  first  principle  theory  to  the  experiments  of 
Hayes  et  al.  in  order  to  assess  the  accuracy  of  their  evalua¬ 
tion  of  the  distribution  function.  We  find  that  overall,  their 
experiments  will  resolve  the  distribution  function  if  reflec¬ 
tions  at  the  collector  barrier  can  be  neglected.  Also,  potential 
fluctuations  due  to  the  discrete  nature  of  impurities  must  be 
negligible  if  their  conclusions  on  the  scattering  rate  are  to  be 
correct. 

In  order  to  precisely  simulate  electron  transport  in  de¬ 
vices  having  heavily  doped  regions,  our  current  model  in¬ 


cludes  the  following  features  (in  addition  to  the  convention¬ 
al  mechanisms ) :  ( i )  elastic-ionized-impurity  scattering,  ( ii ) 
long-range  plasmon  scattering  in  addition  to  short-range 
electron-electron  scattering,  (iii)  coupled  plasmon/phonon 
scattering,  and  (iv)  the  Pauli  exclusion  principle. 

Since  ionized  impurity  scattering  and  the  pair  electron- 
electron  interaction  are  well  known  and  have  been  described 
previously  in  connection  with  Monte  Carlo  methods"  we 
concentrate  here  on  the  addition  of  the  plasmon  effect.  Plas- 
mons  and  coupled  plasmon/phonon  phenomena  have  been 
treated  extensively  in  the  literature.10,13'14  Plasmons  are  due 
to  the  long-range  Coulombic  interaction  between  the  elec¬ 
trons.  They  are  typically  of  long  wavelength  and  therefore 
all  treatments  known  to  the  authors  imply  the  extension  of 
the  electron  gas  over  considerable  distances.  Naturally  in 
ultrasmall  devices  these  long  distances  may  not  be  present. 
In  fact,  a  base  of  600  A  length  doped  with  10l8-cm~3  elec¬ 
trons  will,  on  average,  only  contain  six  electrons  over  its 
length.  These  six  electrons  will  not  support  plasmon  oscilla¬ 
tions  of  wavelength  larger  than  the  base  width.  There  is  also 
another  natural  limit  of  the  plasmon  wavelength  A  as  A  goes 
to  zero. 13  This  limit  is  usually  taken  to  be  twice  the  inverse  of 
the  Thomas-Fermi  wave  vector.13  When  calculating  the 
electron  plasmon/phonon  interaction  we  therefore  include 
these  two  limits,  the  device  length  (as  A— *■»),  and  the 
Thomas-Fermi  wave  vector  (as  A— *0)  in  the  integration  of 
the  scattering  rate.  The  electrons  still  can  interact  in  princi¬ 
ple  with  the  two-dimensional  plasmons  which  propagate 
parallel  to  the  emitter-base  interface.  However,  the  electrons 
entering  the  base  have  initially  negligible  momentum  paral¬ 
lel  to  this  interface.  Also,  as  will  be  described  later,  we  are 
mainly  interested  in  forward  momentum  changes  and  there¬ 
fore  will  ignore  the  coupling  to  the  two-dimensional  pias- 
mons. 


2125 


J.  App).  Phys.  55  (0).  15  March  1986 


0021 -8979/86/0621 25-04502  40 


©  1986  American  Institute  of  Physics 


2125 


New  Ultrafast  Switching  Mechanism  In  Semiconductor  Heterostructures 

K.  Hess,  T.  K.  Higman,  M.  A.  Emanuel  and  J.  J.  Coleman 

Center  for  Compound  Semiconductor  Microelectronics 

Coordinated  Science  Laboratory 

University  of  Illinois 

1406  West  Green  Street 

Urbana,  Illinois  61801 

(Received 

Abstract 

A  new  switching  mechanism  in  a  two-terminal  semiconductor  heterolayer 
structure  is  proposed  which  capitalizes  on  non-linear  electron  temperature 
effects  in  adjacent  heteroiayers.  The  estimated  switching  speed  of  an  optimized 
heterostructure  hot  electron  diode  should  be  extremely  fast,  perhaps  as  fast  as 
200  femtoseconds.  Data  are  presented  on  prototype  devices  which  show  the 
expected  negative  differential  resistance  and  indicate  that  the  basic  physical 
model  is  correct. 
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Intersubband  phonon  overlap  integrals  for  AlGaAs/GaAs  single-well  heterostructures  ,  / 
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The  five  lowest  subband  coupling  coefficients  H^iQ)  are  calculated  for  polar-optical-phonon  scattering 
in  AlojGao^As/Ga As  single-well  heterostructures.  The  wave  functions  are  calculated  self-conststently 
within  the  effective-mass  approximation  including  exchange  and  correlation  effects. 


Much  effort  has  been  made  to  advance  the  understanding 
of  the  electronic  transport  properties  of  q ussi-two- 
dimensional  systems.1  Following  the  pioneering  studies  of 
Stern  and  Howard  for  inversion  layers  of  metal-oxide  semi¬ 
conductor  field-effect  transitors  (MOSFET’s),2  Ando  per¬ 
formed  detailed  studies  of  MOSFET’s  and  AlGaAs/GaAs 
single-well  structures.2  Simultaneously,  lattice  scattering  for 
electrons  in  heterolayers  has  been  extensively  studied.4-4 
Intrasubband  scattering  is  well  understood.  Intersubband 
scattering  has  been  investigated  in  less  detail.4  It  is  the  pur¬ 
pose  of  this  Communication  to  present  an  accurate  calcula¬ 
tion  of  the  overlap  integrals  which  are  necessary  to  deter¬ 
mine  intersubband  scattering  by  phonons.  The  two- 
dimensional  carrier  scattering  rate  between  mth  and  nth 
subbands  is  proportional  tos 

.  (1) 

where  M(Q.q)  is  the  corresponding  three-dimensional  ma¬ 
trix  element,  and  Q  and  q  are  phonon-wave-vector  com¬ 
ponents  parallel  ( xy  plane)  and  normal  (z  direction)  to  the 
layer  interfaces,  respectively.  The  overlap  integral  is  ex¬ 
pressed  as2 

lmn(q  >-  JV„(z  )Fm(z)txp(iqz)dz  ,  (2) 


quantized  energy  levels  ( Em )  but  also  the  corresponding 
wave  functions  [F,(z )  j  for  a  realistic  AlGaAs/GaAs 
single-well  heterostructure.  Using  these  wave  functions,  ac¬ 
curate  values  of  //«.(&)  can  be  obtained.  The  self- 
consistent  calculation  is  briefly  described  in  the  following. 

The  normalized  wave-function  Fm(z)  for  the  mth  sub¬ 
band  can  be  obtained  from  the  Schrodinger  equation 

-—r  )Fm (z )  -  EmFm (z )  .  (4) 

2m  dzl 

The  effective  potential  V(z)  is  given  by  ] 

F(z)--e*(z)+Mz)+Ktc(z)  ,  (5)  ; 

where  ^(z)  is  the  electrostatic  potential  given  by  the  solu¬ 
tion  of  Eq.  (6)  below,  Vk  (z )  is  the  step  function  describing 
the  interface  barrier,  and  V„(z)  is  the  local  exchange  corre¬ 
lation  potential.1  We  take  into  account  the  five  lowest  sub¬ 
bands.  Then,  Poisson's  equation  reads 

- —  I  2  Af/F,2(z)+^(z)-A(0(z)|  ,  (6) 

dz 2  J 

where  N,  represents  the  number  of  electrons  in  subband  / 
and  is  given  by 


where  Fm(z )  is  the  normalized  envelope  function  corre¬ 
sponding  to  the  quantized  energy  of  Em. 

Polar-optical- phonon  scattering  is  one  of  the  dominant 
scattering  mechanisms  above  77  K  in  III- V  compound 
semiconductors.  For  polar-optical-phonon  scattering, 
\M(Q,q  )|2  is  proportional  to  l/((?2  +  q2).  Accordingly,  the 
scattering  rate  is  proportional  to  H„(,Q)/Q,  where  the  mul¬ 
tisubband  coupling  coefficients  are  expressed  as4 

//«.((?)- / /‘ki<fc2Fw,(z,)f*l(zj)exp(-Q|z|-Z2l)  . 

(3) 

Here  Fmiz )  -  Fm(z  )F,(z ). 

Since  most  of  the  electrons  populate  the  lowest  subband 
at  low  temperatures,  the  value  of  H\\(Q)  has  been  worked 
out  in  detail.7  However,  the  electron  occupation  probability 
of  higher  subbands  increases  in  high  electric  fields  (the  elec¬ 
trons  become  hot),  and  both  intrasubband  scattering  in 
higher  suboands  and  intersubband  scattering  among  higher 
subbands  are  no  longer  negligible,  even  at  a  lattice  tempera¬ 
ture  of  77  K.  Under  these  conditions,  the  knowledge  of  the 
coefficients  H^iQ)  is  necessary  in  order  to  accurately  cal¬ 
culate  the  optical  phonon  scattering  rates  for  each  subband. 
Self-consistent  numerical  calculations  provide  not  only  the 


N , 


m*kBT 

irk1 


In 


1  +exp 


Ff~E< 

kgT 


(7) 


Na(z)  and  /VD(z )  are  the  position-dependent  acceptor  and 
donor  concentrations,  and  Ef  is  the  Fermi  energy. 

In  the  calculation,  we  used  an  effective  mass 
m*-0.067mo  and  a  dielectric  constant  e  -  13.0  uniformly  in 
both  GaAs  and  AlGaAs.  We  have  assumed  abrupt  inter¬ 
faces,  a  barrier  height  of  0.3  eV,  and  we  have  neglected  the 
image  force  term  in  Eq.  (S).  All  of  these  approximations 
are  very  reasonable  according  to  the  results  reported  by 
Stem.’  Moreover,  the  mobile  charge  in  the  AlGaAs  layer 
was  assumed  to  be  depleted  except  for  the  penetration  of 
the  wave  function  of  electrons  from  GaAs  into  AlGaAs.  To 
solve  the  Schrodinger  equation,  we  u  3d  the  Numerov 
method,10  which  is  efficient  to  obtain  the  pairs  of  eigen¬ 
values  and  eigenfunctions  accurately  within  short  computa¬ 
tional  times. 

Figure  1  shows  the  calculated  effect  potential  according  to 
Eq.  (5)  together  with  the  five  lowest  subband  energy  levels 
at  77  K.  The  broken  line  indicates  the  electron  concentra¬ 
tion  which  is  derived  from  the  first  term  in  the  right 
parentheses  of  Eq.  (6).  The  structure  parameters  are  shown 
in  the  figure  caption.  The  calculated  energy  differences 


21  6872 


®1985  The  American  Physical  Society 


PHYSICAL  REVIEW  B  VOLUME  33.  NUMBER  8  IS  APRIL  1986 

Monte  Carlo  study  of  electronic  transport  in  AA\ _x  Gax  As/GaAs 
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A  study  of  electronic  transport  in  Al^Ga* As/GaAs  single-well  structures  including  multisub¬ 
band  conduction  at  77  and  300  K  has  been  performed.  The  electronic  states  of  the  quantum  well  are 
calculated  self-consistently  taking  the  five  lowest  subbands  into  account.  The  numerically  obtained 
wave  functions  and  energy  levels  are  used  to  obtain  the  major  two-dimensional  scattering  rates  in 
each  subband.  Polar  optical-  and  acoustic-phonon  (via  deformation-potential)  scattering  are  con¬ 
sidered  including  intersubband  transitions.  For  ionized  impurity  scattering,  the  screening  effects 
due  to  the  five  lowest  subbands  are  taken  into  account  to  obtain  the  Fourier-transformed  Coulomb 
potential.  The  steady-state  and  transient  behavior  of  the  electrons  in  the  well  are  studied  through  a 
Monte  Carlo  particle  simulation.  It  is  shown  that  high  transient  velocities  [(3— 8)X  IQ7  cm/sec]  can 
be  expected  at  low  and  intermediate  fields. 


I.  INTRODUCTION 

Enhanced  electron  mobility  in  modulation-doped 
Alj -jGa, As/GaAs  heterostructures  has  been  suggested 
by  Esaki  and  Tsu1  and  demonstrated  independently  by 
Dingle,  Stormer,  Gossard,  and  Wiegmann.2  Subsequently, 
much  attention  has  been  focused  on  modulation-doped 
heterostructures  in  order  to  realize  ultra-high-speed  field- 
effect  devices,  such  as  the  high-electron-mobility  transis¬ 
tor  (HEMT).3  Among  some  other  combinations  of  III-V 
compound  semiconductors,  the  largest  effort  has  been  de¬ 
voted  to  study  the  Al,_xGax  As/GaAs  system.4 

In  this  system,  the  Al|_xGaxAs  layer  is  doped  n  type, 
while  the  adjacent  GaAs  layer  is  grown  as  pure  as  possible 
(not  intentionally  doped).  The  free  electrons  are 
transferred  to  the  high-purity  GaAs  layer  where  they  pop¬ 
ulate  a  narrow  potential  well  and  form  a  quasi-two- 
dimensional  electron  gas  (Q2D  EG).  Since  the  electrons 
are  spatially  separated  from  the  ionized  donors,  they  ex¬ 
hibit  a  high-electron  mobility  even  at  high  carrier  densi¬ 
ties.  In  fact,  very  impressive  high-mobility  values  have 
been  reported.5-7  It  has  also  been  demonstrated  that  sin¬ 
gle  period  heterostructures  exhibit  better  low-field  mobili¬ 
ties  than  multiple  heterolayers.8 

Not  only  experimental  but  also  theoretical  work  has 
been  stimulated  to  understand  electronic  transport  in 
heterolayers.  Investigations  of  quasi-two-dimensional  sys¬ 
tems9  have  originally  been  performed  for  the  inversion 
layers  of  metal-oxide-semiconductor  field-effect  transis¬ 
tors  (MOSFET’s).  Following  the  pioneering  studies  of 
Stem  and  Howard  for  MOSFET’s,10  Ando  carried  out  de¬ 
tailed  calculations  of  MOSFET’s  and  Al|_xGax  As/GaAs 
single-well  structures."  In  these  studies,  the  subband 
structure  has  been  calculated  using  wave  functions  and 
energy  levels  which  are  obtained  by  a  self-consistent  cal¬ 
culation  based  upon  the  effective-mass  approximation. 

Simultaneously,  lattice  scattering12- 16  and  impurity 
scattering14,16-19  for  electrons  in  heterolayers  have  been 
extensively  investigated.  Many  of  the  above  treatments 
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have  assumed  square-well  potentials  and  approximate 
wave  functions  to  express  the  scattering  rates  analytically. 
Walukiewicz,  Ruda,  Lagowski,  and  Gatos20  used  a 
triangular-well  approximation  to  discuss  the  electron  mo¬ 
bility  more  realistically.  For  impurity  scattering  in  Si  in¬ 
version  layers,  Mori  and  Ando  reported  formulas  using 
variational  functions.17  However,  they  focused  their  at¬ 
tention  upon  low  temperatures  to  discuss  the  ultimate  lim¬ 
its  of  the  mobility  values.  Their  final  results  are  not  appl¬ 
icable  at  high  temperatures. 

In  all  previous  publications,  the  electron  mobility  values 
have  been  calculated  from  a  relaxation  time  which  is 
closely  related  to  the  inverse  of  the  scattering  rate.  For 
device  applications,  we  need  to  understand  the  warm- 
and/or  hot-electron  behavior  of  the  quantum  well  system 
at  and  above  77  K.  To  study  the  transport  properties  in 
this  regime,  the  Boltzmann  transport  equation  must  be 
solved  taking  into  account  the  scattering  in  each  subband 
of  the  Al|_j,GaxAs/GaAs  single-well  potential.  A 
Monte  Carlo  method21  can  bypass  the  difficulties  in 
directly  solving  the  complicated  system  of  equations. 
Hot-electron  velocity  characteristics  for  a  strictly  two- 
dimensional  system  confined  in  a  square  well  of 
Al|_xGax  As/GaAs  have  already  been  studied  previously 
using  the  Monte  Carlo  method.22 

In  this  paper,  we  present  more  precise  calculations  of 
electron  transport  properties  for  Al,_xGax As/GaAs 
single-well  heterostructures  including  multisubband  con¬ 
duction  by  means  of  an  ensemble  Monte  Carlo  method. 
In  Sec.  II,  the  electronic  states  of  the  quantum  well  are 
calculated  self-consistently  taking  the  five  lowest  sub¬ 
bands  into  account.  The  numerically  obtained  wave  func¬ 
tions  and  energy  levels  are  used  to  calculate  the  major 
two-dimensional  scattering  rates  at  77  and  300  K  in  Sec. 
III.  Polar  optical- phonon  scattering  and  acoustic-phonon 
scattering  (via  the  deformation  potential)  are  studied  for 
each  subband  including  intersubband  transitions.  For  ion¬ 
ized  impurity  scattering,  the  screening  effects  due  to  the 
five  lowest  subbands  are  considered  to  calculate  the 
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A  time-dependent  ensemble  Monte  Carlo  simulation  has  been  developed  which  includes  the 
computation  of  the  self-consistent  electric  field  and  electron  concentration.  The  Monte  Carlo 
code  is  used  to  study  transient  and  steady-state  transport  in  planar-doped-barrier  (PDB)  GaAs 
diodes  and  transistors.  Results  are  presented  for  the  temporal  and  spatial  evolution  of  electronic 
transport  through  the  PDB.  We  show  that  (i)  average  velocities  achievable  with  a  PDB  are 
significantly  above  the  steady-state  velocities  by  about  a  factor  of  3,  (ii)  the  rise  time  of  the  electric 
current  in  the  PDB  is  in  the  subpicosecond  range  for  diode  length  below  2000  A.,  (iii)  a  pronounced 
displacement  current  dominates  the  transient  electronics  at  onset  of  injection,  and  (iv)  polar 
optical  phonon  absorption  contributes  significantly  to  the  thermionic-emission  current  over  the 
barriers. 


I.  INTRODUCTION 

Ultrathin  GaAs  planar-doped  structures  have  attracted 
considerable  attention  in  recent  years.'  They  provide  an  al¬ 
ternative  to  Schottky  barriers  due  to  their  fast  switching, 
mixing,  and  rectifying  properties  and  are  useful  for  hot-elec- 
tron  injection  in  connection  with  nonstationary  transport  in 
submicron  devices.  In  the  past,  several  theoretical  analyses 
have  been  performed  ranging  from  highly  simplified  closed- 
form  analytical  treatments1'3  to  complicated  Monte  Carlo 
simulations.6  The  earliest  analytical  treatment  was  present¬ 
ed  by  Malik,'  who  made  use  of  the  Richardson  thermionic- 
emission  theory.  Kazarinov  and  Luryi  analyzed  the  device 
operation  by  using  the  drift-diffusion  model2;  Cook  replaced 
the  assumption  of  a  constant  mobility  with  an  energy-trans¬ 
port  model4  based  on  Stratton’s  formalism7;  and  Littlejohn 
et  al.  presented  a  Monte  Carlo  model  with  fixed  potential  to 
study  velocity  overshoot.6  These  studies  have  demonstrated 
that  hot-electron  effects,  transient  phenomena,  and  the  self- 
consistency  of  the  electric  field  are  crucial  in  the  operation  of 
planar-doped-barrier  devices.  In  this  paper,  we  therefore 
have  developed  a  Monte  Carlo  code  to  include  all  these  fea¬ 
tures;  in  addition,  we  have  included  in  our  simulation  the 
interactions  of  the  electrons  among  themselves,  although 
this  interaction  is  of  little  significance  for  the  structures  ana¬ 
lyzed  in  this  paper.  Using  our  Monte  Carlo  simulations,  we 
demonstrate  and  analyze  the  significance  of  velocity  over¬ 
shoot  and  propose  a  novel  structure  which  extends  velocity 
overshoot  over  a  submicron  spatial  scale. 

II.  METHOD  OF  ANALYSIS 
A.  D«vic«  structure 

A  planar-doped  structure  consists  of  single  or  multiple 
layers  which  are  highly  doped  and  sandwiched  between  in¬ 
trinsic  material.  By  controlling  doping  concentrations  and 
positions  of  these  layers,  one  can  tailor  the  conduction  band 
profile.  Typical  examples  of  planar-doped  devices  are  n*-i- 
p+-i-n  +  diodes  and  hot-electron  transistors.8  A  schematic 


band  diagram  for  an  n +-i-p+-i-n  +  diode  is  shown  in  Fig.  1. 
Notice  that  the  sandwiched  p+  layer  is  thin  enough  to  be 
totally  depleted  in  equilibrium. 

B.  Monte  Carlo  model  and  boundary  conditions 

Our  ensemble  Monte  Carlo  method  incorporates  the  A 
L-X  valleys  of  the  GaAs  conduction  band  and  includes  ef¬ 
fects  of  nonparabolicity.  Materia]  parameters  such  as  non- 
parabolicities,  valley  separations,  and  transverse  and 
longitudinal  effective  masses  are  obtained  from  experimen¬ 
tal  results  as  available  and  from  pseudopotential  calcula¬ 
tions.9  The  scattering  mechanisms  included  are  electron- 


Position  (A) 

FIG.  I .  Schematic  representation  of  the  band  diagram  of  the  n  *  - i-p  *  -i-n  * 
structure  shown  in  the  inset  of  the  figure  with  a  bias  voltage  of  0.82  V.  The 
area  doping  concentration  in  the p *  plane  is  l.Ox  1011  cm  ~ !  and  the  vol¬ 
ume  concentration  in  the  n *  contact  region  is  l.Ox  10"  cm “ ’ 
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Electron  transport  in  GaAs/AIGaAs  superlattices  was  investigated  by  means  of 
Monte  Carlo  simulations  The  sources  of  scattering  considered  are  acoustic  and  polar- 
optical  phonons  and  ionized  impurities  Collision  broadening  of  the  extended  slates  is 
taken  into  account  and  ns  effect  on  the  transport  parameters  is  assessed. 


1.  Introduction 

The  concept  of  semiconductor  superlallices  has  been 
introduced  a  lew  years  ago  by  Hsaki  and  'Hsu  [l]  with  the 
expectation  that  Bloch  oscillations,  a  rather  prohibitive 
phenomenon  in  bulk  material,  may  really  occur  in  such 
artificial  structures  In  superlallices.  the  artificial  periodi¬ 
city  creates  "mini"  Bnllouin  zones  with  boundaries  that  can 
be  reached  by  the  electron  betore  scattering  to  some  other 
part  of  the  zone. 

The  shape  of  the  energy  dispersion  relation  [  1  ]  shows 
that  as  the  electron  has  the  chance  to  cross  the  minizone  it 
will  reach  a  region  where  ns  mass  becomes  negative  allow¬ 
ing  for  the  observation  of  negative  differential  resistance. 

The  existence  of  Bloch  oscillations  and  negative 
differential  resistance  has  been  predicted  using  the  assump¬ 
tion  of  a  constant  scattering  time  [1.2]  and  also  by  Monte 
Cario  calculations  [3.4]  of  the  drill  velocity  when  an  elec¬ 
tric  field  parallel  to  the  growth  axis  of  the  superlallice  is 
applied.  These  calculations  have  been  performed  by  assum¬ 
ing  that  the  electrons  populate  only  the  first  superlallice 
subband  and  that  the  scattering  agents  are  acoustic  and 
polar  optical  phonons. 

It  is  the  purpose  of  this  paper  to  investigate  the 
influence  ol  collision  broadening  on  superlallice  transport. 
The  reason  lor  our  investigation  is  the  small  energy-width 
ol  the  superlallice  subbands  (<  50meV).  We  therefore 
simulated  steady  stale  transport  in  a  superlallice  by  going 
beyond  the  semictassical  Boltzmann  equation  along  lines 
previously  used  by  Chang  el  al  [6], 

We  also  include  as  an  additional  scattering  mechanism 
of  our  simulation  ionized  impurity  scattering. 


II.  Model  of  the  Superlallice  Studied 

Our  model  is  an  idealization  of  a  GaAs-Al  Ga.  As 
superlallice  having  a  mole  fraction  of  x  ^  0.3.  Given  low 
x-values  and  the  rather  low  values  of  the  applied  field  it 
was  assumed,  as  done  by  previous  aulhorsJ'',  that  the  elec¬ 


trons  will  not  transfer  to  the  subsidiary  L  and  X  minimum 
of  GaAs  or  AlxGaj_xAs. 

The  energy  dispersion  relation  can  therefore  be 
approximated  by: 

£( k  )=£„+i(-l)"  A£„cos*,rf +-5l_(*,J+*,2)  (i) 

4  2m 

where  the  superlattice  growth  direction  is  the  z-axis.  (£„ 
is  the  center  energy  of  the  nth  subband  and  &£„  its  width.) 

in  determining  the  values  of  E„  and  .  we  used  the 
approach  of  Leburton  and  Kahen  [7],  Their  method  is  an 
extension  of  a  simple  Kronig-Penney  periodic  potential  to 
account  lor  some  details  of  a  more  realistic  electron  struc¬ 
ture.  Using  their  results  for  a  superlallice  with  x  »  0.1  ./.„ 
(well  width)  *  LH  (barrier  width)  *  8a„  (alP  being  the  GaAs 
lattice  constant )  and  assuming  a  0.65  0.35  rule  for  the  band 
edge  discontinuity  ol  the  conduction  and  valence  bands,  we 
obtain: 

£,=  0.0721  eV  A£,=  0.01726  eV 

£\=  0.2109  eV  A£,=  0.0551  eV 

The  lorbtdden  zone  between  the  first  and  second  subband  is 
then  0.1289  eV. 


ill.  Simulation  method 

The  transport  parameters  were  calculated  with  the 
standard  one  electron  steady-slate  Monte  Carlo  simulation 
whose  details  have  been  *  described  in  detail  by  many 
authors  (see  for  example  Fawcett  et  al.  [5]).  The  electron 
was  launched  at  the  bottom  of  the  miniband  and  its  trajec¬ 
tory  monitored  m  the  first  superlallice  minizone.  Flight 
times  and  the  types  ol  scattering  mechanisms  ending  the 
electron  drill  were  determined  stochastically  according  to 
the  scattering  probabilities. 
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Transient  transport  of  holes  in  GaAs  is  studied  under  the  conditions  of  high  energy  injection.  This 
study  is  made  using  a  Monte  Carlo  simulation  with  the  unique  inclusion  of  a  realistic  band 
structure  based  on  the  K  *P  method.  The  results  reported  herein  show  that  a  significant  velocity 
overshoot  occurs  at  low  applied  electric  fields  over  distances  greater  than  1 500  A.  As  the  applied 
field  increases,  the  effect  of  the  overshoot  upon  the  transit  time  becomes  less  pronounced.  There  is 
no  disceraable  gain  in  the  drift  velocity  for  holes  injected  at  high  energy,  0.1  eV,  as  compared  to 
holes  injected  at  zero  launching  energy  for  the  distances  considered  here. 

PACS  numbers:  72.20.Js 


Much  attention  has  been  focused  recently  upon  deve¬ 
loping  high  speed  devices  based  upon  velocity  overshoot.1-3 
As  originally  proposed  by  Ruch,6  use  of  velocity  overshoot 
of  electrons  in  GaAs  and  other  III-V  compounds  can  give  a 
significant  increase  in  carrier  velocity  over  the  steady-state 
value.  The  recent  work  of  Tang  and  Hess7  and  Brennan  et 
al.*  has  shown  that  there  exists  a  limited  set  of  conditions, 
the  collision  free  window,  under  which  velocity  overshoot 
can  be  appreciable  for  electrons  in  GaAs  and  InP.  However, 
the  distances  over  which  velocity  overshoot  persists  have 
been  found  to  be  small.7,8  There  may  be  ways  of  artificially 
lengthening  the  range  over  which  velocity  overshoot  can  be 
attained.  One  method  suggests  that  repeated  overshoot  can 
be  achieved  with  electrons  in  GaAs  through  the  use  of  poten¬ 
tial  steps. 9,10  The  recent  work  of  Brennan  and  Hess 10  shows 
that  velocity  overshoot  of  electrons  can  be  maintained  over 
distances  of  0.5  ft  or  more  by  use  of  staircase  structures. 

Ruch  found  that  the  electron  velocity  overshoots  its  sat¬ 
uration  value  in  both  silicon  and  gallium  arsenide.6  How¬ 
ever,  the  velocity  overshoot  persists  an  order  of  magnitude 
longer  in  GaAs  than  in  Si  due  in  part  to  the  much  faster 
energy  relaxation  in  Si  at  low  energies.  The  electron  drift 
velocity  also  is  substantially  higher  in  GaAs  than  in  Si  at  low 
applied  fields1 1,12  because  of  the  different  effective  masses  of 
the  elecrons  in  each  material.  In  GaAs,  polar  optical  scatter¬ 
ing  is  the  dominant  energy  relaxation  mechanism  which  on 
average  is  weaker  than  the  deformation  potential  scattering 
present  in  Si.  Therefore,  it  takes  longer  for  the  electron  over¬ 
shoot  to  relax  in  GaAs  than  in  Si,  leading  to  much  higher 
electron  drift  velocities  over  longer  distances.7 

The  presence  of  heavy  holes  and  deformation  potential 
scattering  in  the  valence  band  of  GaAs  at  low  energy  greatly 
diminishes  the  hole  velocities.  Therefore,  one  would  expect 
that  velocity  overshoot  of  holes  would  be  limited  by  the 
strong  energy  relaxation  much  as  it  is  in  silicon.  However, 
our  calculations  show  that  the  velocity  overshoot  of  holes  in 
GaAs  is  significant  under  certain  conditions  of  field  and 
launching  energy.  Velocity  overshoot  of  holes  in  GaAs  is  not 
as  pronounced  as  that  of  electrons  in  GaAs  but  it  is  larger 
than  for  electrons  in  silicon,  as  calculated  by  Ruch.6 


Our  calculations  are  performed  using  a  Monte  Carlo 
program  with  a  realistic  valence  band  structure13  based  on  a 
K  *P  calculation. 14  It  is  essential  that  a  realistic  band  struc¬ 
ture  be  used  because  the  valence  bands  are  strongly  warped 
and  nonparabolic  even  at  low  energy. 13  We  include  the  split- 
off  band  as  well  as  the  heavy  and  light  hole  bands  in  the 
calculation.  The  hole-phonon  scattering  rate  is  determined 
using  a  field  theoretic  approach16,17  and  is  based  on  the  total 
density  of  states  present  in  all  the  bands.  The  probability  of  a 
hole  transfering  from  one  band  to  another  via  interband  scat¬ 
tering  is  determined  by  directly  calculating  the  overlap  inte¬ 
grals.13 

The  predominant  scattering  mechanisms  are  polar  opti¬ 
cal  and  deformation  potential  scattering. 18,19  The  principal 
scattering  agents  are  acoustic,  nonpolar,  and  polar  optical 
phonons.  The  acoustic  scattering  rate  is  calculated  using  the 
method  of  Canali  et  al.,20  and  the  value  of  the  acoustic 
phonon  coupling  constant  E  f  is  determined  from  the  work 
of  Costato  and  Reggiani. 18  The  value  of  the  deformation  po¬ 
tential  constants  “a,  b,  d,"  which  are  used  in  determining 
£ 2  ,  are  taken  from  Wiley.21  The  transition  probabilities  for 
nonpolar  optical  scattering  can  be  written  in  a  manner  that  is 
completely  analogous  to  that  for  acoustic  phonon  scatter¬ 
ing22  containing  the  optical  deformation  potential  constant 
However,  it  is  useful  to  express  the  nonpolar  optical  scat¬ 
tering  rate  in  terms  of  the  optical  phonon  coupling  constant 
[DK  )2.  (DK  )2  can  easily  be  determined  from  £?.21  All  of  the 
parameters  used  in  determining  the  hole-phonon  scattering 
rate  will  be  collected  in  a  future  publication.23 

The  overall  scattering  rate  is  calculated  using  both  the 
results  of  Costato  and  Reggiani18,19  and  a  field  theoretic  ap¬ 
proach.16,17  At  low  energy,  the  total  phonon  scattering  rate 
is  calculated  from  the  sum  of  the  acoustic,  nonpolar,  and 
polar  optical  scattering  rates  as  determined  by  Costato  and 
Reggiani.18,19  Above  approximately  an  energy  of  0.5  eV,  the 
total  scattering  rate  is  calculated  based  on  the  total  density  of 
states  of  all  three  valence  bands.  This  rate  is  adjusted  to 
match  the  low  energy  scattering  rate  at  0.5  eV.  In  the  calcu¬ 
lations  presented  here,  the  self-energy  effects  are  significant 
at  high  applied  fields:  100  and  200  kV/cm. 
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The  total  scattering  rate  and  the  transition  probability  for  electron-phonon  interaction  in  1-D  and 
2-D  semiconductor  materials  are  calculated  in  taking  into  account  the  finite  dimensions  of  the 
structure.  Although  noticeable,  size  effects  on  the  scattering  rate  are  generally  small,  with  more 
pronounced  features  for  1-D  structures  than  for  2-D  structures.  For  2-D  layers,  our  theory  agrees 
with  recent  experimental  results  whereas  it  contradicts  the  previous  theory  predicting  large  sue 
effects  and  mass-independent  electron-phonon  scattering  rates.  In  1-D  structures  singularities  in 
the  phonon  emission  rate  appear  as  a  natural  consequence  of  the  1-D  density  of  states.  However, 
for  high  energy  the  1-D  emission  rate  is  found  smaller  than  the  corresponding  3-D  rate.  An 
additional  consequence  of  the  confinement  is  the  quenching  of  the  phonon  absorption  rate. 


I.  INTRODUCTION 

The  remarkable  developments  of  fine  line  lithography 
and  the  new  epitaxial  technologies  of  molecular  beam  epi¬ 
taxy  (MBE)  and  metalorganic  chemical  vapor  deposition 
(MOCVD)  have  added  astonishing  possibilities  to  the  fabri¬ 
cation  of  new  artificial  materials.  Alternating  structures  of 
lattice  matched  III-V  semiconductor  compounds,  as  for  ex¬ 
ample  AlxGa,  _,As  are  presently  fabricated  with  typical 
thicknesses  of  10  nm  or  less  and  with  accurately  controlled 
doping.  In  these  artificial  materials,  the  layer  dimensions  are 
comparable  to  the  de  Broglie  wavelength  A,  and  the  elec¬ 
trons  are  confined  in  the  ultrafine  structure.  Transport  is 
essentially  two  dimensional  (2-D)  and  presents  new  charac¬ 
teristics  not  exhibited  in  bulk  semiconductors.  For  example, 
the  mobility  parallel  to  the  plane  of  the  layer  has  been  shown 
to  be  improved  drastically  by  modulation  doping  of  the 
structures  and  values  as  high  as  106  cmWs  have  been 
achieved. 1  Recently,  an  alternative  structure  based  on  the 
confinement  of  electrons  in  a  “wire”  semiconductor  has 
been  proposed  to  improve  further  the  conducting  properties 
of  the  materials.2  In  these  one-dimensional  ( 1-D)  III-V  syn¬ 
thetic  materials,  the  electron  gas  is  quantized  in  two  trans¬ 
verse  directions  and  the  charge  carriers  can  only  move  in  the 
longitudinal  direction.  GaAs-AlA?  quantum  well  wires  have 
already  been  realized  by  Petroff  et  al.1  and  other  structures 
based  on  field  effects  promise  to  be  feasible  in  the  near  fu¬ 
ture.4,5 

Because  of  the  limited  number  of  available  final  states 
during  the  scattering  process  (forward  or  backward  scatter¬ 
ing),  the  1-D  electron  mobility  is  seemingly  enhanced  and  its 
value  has  been  estimated  to  be  far  beyond  106  cmW  s.6  The 
obvious  potential  of  these  new  structures  for  high-speed  de¬ 
vices  makes  the  knowledge  of  their  transport  parameters  de¬ 
sirable.  Although  transport  theory  in  1-D  semiconductor 
structures  is  still  in  its  infancy,  calculations  of  the  most  com¬ 
mon  scattering  mechanisms  impurity,6-7  acoustical 
phonon,*  and  optical  phonon9  have  already  been  performed. 
However,  interaction  with  polar  optical  phonons  (POP)  mer¬ 
its  special  attention;  the  large  energy  exchange  between  car¬ 


riers  and  lattice  vibrations  influences  considerably  transport 
at  room  temperature.  Moreover,  the  POP  interaction  is  a 
major  scattering  mechanism  responsible  for  hot  electron  ef¬ 
fects.  >0  So  far,  the  dependence  of  the  POP  scattering  on  the 
“wire”  cross  section  has  not  been  investigated.  Yet,  wire  pa¬ 
rameters  appear  to  be  technologically  controllable,  and 
owing  to  device  applications,  it  is  of  importance  to  appre¬ 
ciate  realistically  the  confinement  effects  on  the  electron- 
POP  interaction.  In  this  paper  we  derive  the  total  POP  scat¬ 
tering  rate  in  1-D  semiconductor  structures  and  discuss  the 
influence  of  size  effects  on  the  electron-phonon  interaction. 
Furthermore,  we  reconsider  the  corresponding  2-D  transi¬ 
tion  probability.  2-D  electron-POP  interaction  represents  an 
important  limitation  on  the  high-speed  performance  of  High 
Electron  Mobility  Transistor  (HEMT).  Although  experi¬ 
mental  data  show  a  weakening  of  the  electron-phonon  inter¬ 
action  in  2-D  structures  compared  to  3-D  materials"12 
theoretical  models  predict  a  larger  scattering  rate. 13  This  dis¬ 
crepancy  has  been  attributed  to  the  screening  of  the  electron- 
POP  interaction  by  high-density  carriers  in  2-D  structures. 
However,  we  believe  that  the  discrepancy  is  partially  due  to 
crude  approximations  previously  used  for  the  2-D  calcula¬ 
tion,  e.g.,  (strictly  2-D  gas),  where  q,  is  the  trans¬ 

verse  component  of  the  phonon  wave  vector  and  L,  is  the 
layer  width.  '*• 13  Hence,  Ridley,  introducing  the  momentum- 
conservation  approximation  (MCA)  has  obtained  an  expres¬ 
sion  for  the  2-D  scattering  rate  which  is  mass  independent 
and  varies  as  1  /L, . 16  Prior  to  these  results.  Price  obtained  a 
■[in  law  for  the  2-D  collision  rate,  the  same  as  for  3-D  elec¬ 
tron  gases. 17  He  also  has  shown  that  the  size  effects  are 
weaker  than  predicted  by  Ridley.  In  the  following,  we  clarify 
this  question  by  calculating  exactly  the  phonon  scattering 
rate  for  2-D  electron  gases.  Our  theory  generally  agrees  with 
experimental  data  from  Chiu  et  a/.12  and  shows  that  size 
effects  lower  the  POP  scattering  rate.  Screening  effects,  how¬ 
ever,  are  not  considered  here  and  likely  account  for  further 
lowering  of  the  electron-POP  interaction  as  claimed  by  Shah 
et  al."  For  1-D  electron  gas,  size  effects  are  shown  to  have 
important  consequences  on  the  POP  emission  rate. 
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Summary:  The  principles  of  real  space  transfer  and  its  phenomenological  description  in  terms 
of  the  concepts  of  electron  temperature  and  quasi  Fermi  levels  are  reviewed.  It  is  shown  that 
real  space  transfer  is  a  mechanism  to  achieve  ultrafast  switching  and  storage  of  charge  carriers. 
The  real  space  transfer-glow  cat  hod  analogy,  which  demonstrates  the  existence  of  a  new  transis¬ 
tor  principle,  is  discussed  in  detail. 


Introduction 

The  technologies  of  molecular  beam  epitaxy  (MBE)  and  metalorganic  chemical 
vapor  deposition  (MOCVD)  have  opened  a  world  of  possibilities  to  create  new 
forms  of  II1-V  semiconductor  heterolayer  materials.  Structures  hundred  times 
as  sophisticated  and  small  as  the  smallest  feature  sizes  in  current  silicon  technology 
can  be  fabricated  with  ease.  This  fine  tuned  variability  of  structure  and  boundary 
conditions  offers  many  opportunities  for  the  development  of  new  device  concepts 
and  research  in  this  area  is  rapidly  expanding  [  1 ). 

Advantages  of  I1I-V  compounds  with  respect  to  steady  state  mobility  and  transient 
electronic  transport  [1]  (velocity  overshoot)  have  long  been  known.  However, 
the  realization  of  these  advantages  in  applications  have  been  impeded  by  technolog¬ 
ical  difficulties.  MBE  and  MOCVD  have  provided  a  quantum  leap  in  material 
quality  and  controlability,  which  together  with  the  ideas  of  selective  (modulation) 
doping  (2]  and  achievable  abruptness  of  heterointerfaces  have  increased  the  mobili¬ 
ty  values  and  potential  device  speeds  by  an  order  of  magnitude.  Studies  of  modula¬ 
tion  doping  have  resulted  in  successful  designs  of  the  high  electron  mobility  transis¬ 
tor  [3]  which  draws  its  major  advantage  from  the  low  source  access  resistance. 
Simultaneously  to  these  developments  several  other  new  effects,  peculiar  to  the 
boundary  conditions  of  complicated  heterolayer  structures,  have  been  discovered. 
Examples  are  quantum  well  lasers  (4],  and  real  space  transfer  devices  (5, 6].  It  is  the 
purpose  of  this  paper  to  describe  the  physical  basis  of  real  space  transfer  and  iu 
xxxxx  xxxxx 


